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Pionic excitations are analyzed in isospin-asymmetric Fermi gas of constituent 



cn 

. quarks. It is found that in addition to the usual pionic excitations, there exists 

0^ I a very low collective mode with quantum numbers of the charged pion (spin-isospin 

sound). In the chiral limit the excitation energy of one of the modes scales as the 
^ . current quark mass, and the mode saturates current-algebraic sum rules. This is in 

agreement with predictions of current algebra for isospin-asymmetric medium. 



C ■ Recently chiral current algebra has been applied [1,2] to systems with finite 

chemical potentials in order to study excitations in such systems. A formal pre- 
^ I diction made by Cohen and one of us (WB) is that in the isospin-asymmetric 

I nuclear medium there exists an excitation with quantum numbers of the 

charged pion (vr^ for medium with negative isospin density), which in the 
chiral limit becomes very soft [1]. Its excitation energy is proportional to the 
current quark mass m, and not ^/m as in the case of the vacuum, and in this 
sense the mode is softer than the usual pseudo-Goldstone boson.Q If such a 
soft mode exists at the physical value of m in dense nuclear medium, and if it 
is strongly coupled, it would bring important phenomenological consequences 
by dominating excitation functions and thermal properties of the system at 
low energies. Unfortunately, current algebra gives us no hint as to whether 
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the physical value of m is low enough such that the formal results of Ref. [1] 
have practical importance. After all, in a nuclear system we have additional 
energy scales, e.g. chemical potentials or energies of singic-particle excitations, 
and m, or rather the vacuum value of the pion mass, m^r, need not be small 
compared to these scales. 

In order to gain some quantitative insight we study a model in which the asym- 
metric nuclear medium is composed of the Fermi gas of constituent quarks of 
the Nambu-Jona-Lasinio model [3]. Our results are as follows: Firstly, since 
the model obeys all requirements of current algebra, the chiral soft mode 
emerges when m is decreased. However, this mode is not necessarily a contin- 
uation of the vacuum pion branch. Another branch of low excitations exists 
(~ 10 MeV), which is analogous to the spin-isospin sound [4,5] found in con- 
ventional approaches to neutron matter. We show that this mode contributes 
sizably to current- algebraic sum rules, also at the physical value of m. 

Let us begin with a reminder of the sum rules of Ref. [1,2]: 
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where J" q = g7o75|T"g is the time component of the axial vector current, \C) 
denotes a uniform medium subject to constraints imposed externally, | j^) 
are excited states with three-momentum equal zero of isospin ±1, and Ej± are 
their excitation energies. The constraints are the baryon density, pb, and the 
isospin density p/=i. In Eq. (1) we recognize the generalization of the Gell- 
Mann-Oakes-Renner sum rule for finite densities. The "isovector" sum rule 
(2) is nontrivial only for systems with isospin asymmetry. As shown in Ref. [1], 
since the isovector chemical potential remains finite as m — > 0, in a medium 
with negative isospin density only positive-isospin modes can contribute to 
the above sum rules in the chiral limit. This implies, under some additional 
weak assumptions, the existence of a "chiral soft mode" | -|-) for which 
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This mode saturates both sum rules in the chiral limit. On the other hand, 
there are no soft modes with negative isospin, and Ej- ~ 1. The presented 
behavior is radically different from the case of the vacuum or symmetric mat- 
ter, where the lowest 7r+ and 7r~ excitations scale as y/m. The ttq excitations 
scale in all cases as ^/m [1] . 

Now we pass to an illustration of these general results in a simple model: 
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the Nambu-Jona-Lasinio model with scalar and vector interactions [3,6]. The 
model is consistent with current algebra [2], and its variants have been applied 
extensively to describe both meson and baryon physics. The Lagrangian is 

C = q{\ ^-m)q-r^ {{qqf + {qil.r'^qf) + ^ {{qr^qf + {qil.qf) 

- ^ ((?7m^"?)' + (?757.r'^?)') - %(?7.?)' ■ (4) 

Using the Hartree approximation one arrives at self-consistency equations for 
the values of the scalar-isoscalar field the scalar-isovector field 5, the time 
component of the neutral vector-isovector field p, and the time component of 
the vector-isoscalar field uj : 



S = m ~ Gfjiuu + dd), 5 = —Gs{uu — dd), 

p=2Gp{u+u - d+d), uj = G^{u^u + d+d). (5) 

The scalar and vector densities of the u quark are equal to (analogously for 
the d quark) 

/d^k M r d^k 

We have defined 6^ = Q{ku — \k\) — Q{A — \k\), where A is the sharp three- 
momentum cut-off used in this paper, and ku and kd are the u and d quark 
Fermi momenta. We have also introduced scalar self-energies M„ — S + 5 and 
Md — S — S. Self-consistency requires that the quark propagators are evaluated 
with mean- fields (5): 

^uld =i> - ^o(^f + - M^id -f is sgn(/i„/d - po) , (7) 

where /i„ and ^i^ are the chemical potentials of the u and d quarks. Next, we 
consider the charged pion propagator. We choose kd > ku, as for instance in 
neutron matter. Were there no vector-isovector meson effects {i.e. Gp — 0), 
the inverse pion propagator would have a simple form l — G^JpF, where Jpp — 
— iTr / -^ii^jbSulp + \q)^5Sd{p — \q). However, in presence of vector-isovector 
coupling there occurs tt — Ai mixing, and to find excitation energies with pion 
quantum numbers one has to find zeros of the determinant of the inverse tt—Ai 
propagator matrix (see e.g. Ref. [2] for details). In the considered model this 
determinant can be written as 

?o - P I -5 V qoiqo- P) J I Si} 
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where qo is the frequency of the excitation, and 




\i (2^)' (P - Qo)' + 2(p - qo)^k^ + M2 + 4SS 



+ {u^ d, -p, S ^ -6, qo -qo) • 



(9) 



Let us look at the analytic structure of D, or Jpp, in go- The state \C) consists 
of the Fermi seas of d and u quarks, with kd > ku, as well as of the Dirac sea 
occupied down to the cut-off A. A positive-charge Fermi sea excitation moves a 
quark from the occupied d level to an unoccupied u level. Pauli blocking allows 



this when p+^kj + M^ - ^kj + Mj < qo < p + ^kl + - ^kl + M^. 



Thus, within these boundaries D{qo) possesses a cut. There are also additional 
cuts corresponding to quark- ant iquark creation, but their effects are negligible 
in our study. As explained e.g. in [4,5] in the framework of conventional nuclear 
physics, it is possible for the pion propagator in neutron matter to have an 
additional pole at very low excitation energies: the spin-isospin sound. We will 
show that an analogous phenomenon occurs in our model. 

For the numerical study we use the following parameter sets: 

1:G^ = 7.55GeV-^ Gs = 5.41GeV-^ Gp = 7.09GeV-^ 

A = 750MeV, m = 3.52MeV, 
ll :Ga = 4.35GeV"^ Gs = 3.34GeV"^ Gp = 12.4GeV"^ 
A = 954MeV, m = 2.03MeV. 

Set I is the SU(3) fit of Ref. [7], which fits the = 139MeV, = 93MeV, 
rrirj = 519MeV, and rup = 765MeV. This leaves one parameter out of the 
original five undetermined, which is expressed through the constituent quark 
mass in the vacuum, set arbitrarily to M = 361MeV. Set II also fits m7r, F^^, 
rrir^, and M = 361MeV, but not nip. Instead, we chose a much larger value of 
Gp. Such greater values are needed if, for instance, one wishes to fit the an 
pion-pion scattering length [6]. Since the excitations carry no baryon number, 
the value of G^, is irrelevant. 

We define x = ps/po, where pB is the baryon density and po = 0.17fm^'^ is the 
nuclear saturation density. The relative concentration of d quarks is denoted 
by 2/ = Pd/{Pu + Pd)- For the results presented below we set y = 2/3, which 
corresponds to pure neutron matter. With y fixed, we vary the x variable, 
increasing the density from to a few times the saturation density. Figure 
1 shows the results for the parameter set I. In (a) we plot the mean fields 
(5). In (b) we show the excitation energies of what we call the "big" modes, 
the ones that at a; = (vacuum) become the usual pions. The excitation 
energy of n'^ drops to about lOOMeV a,t x — 2, and then increases again at 
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higher x. The excitation energy of vr" increases sharply, which is a result of 
Pauli blocking. At low densities the behavior of these excitation energies is 
controlled by the Weinberg- Tomosawa term in the tt — N scattering [2] . The 
lower part of Fig. 1(b) shows the excitation energy of the third solution of 
the equation D{qo) — 0. This mode is the spin-isospin sound, Trf. Its energy 
is small, rising up to only 6MeV at a; = 2. The shaded area shows the cut 
discussed after Eq. (9). The spin-isospin sound emerges from the cut at a very 
low but finite value of ,t. A similar phenomenon has also been obtained in 
studies of the kaon propagator in nuclear matter [8,9]. How relevant is this 
mode? A natural measure of its strength is provided by the contributions to 
the sum rules (1-2). Explicitly, we have 
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The contributions of the three modes to the sum rules are plotted in Fig. l(c- 
d). Whereas in the GMOR sum rule the spin-isospin sound contributes up to 
10% (which is still noticeable), in the isovector sum rule it dominates over 
the "big" modes at x > 2, and Bit x — 3 practically saturates the sum rule. 
In that sense the mode is strong. Note that the three modes satisfy the sum 
rules at the 99% level, which means that the contributions from cuts are tiny. 
Next, we present a formal study of the chiral limit. In Fig. 2(a) we plot the 
excitation energies of the three modes at a; = 2 and y = 2/3 as functions of 
a, which is the ratio of the current quark mass to its physical value from the 
parameter set I. Hence a — 1 corresponds to the physical case, and a = is 
the strict chiral limit. The spin-isospin sound becomes the chiral soft mode of 
Eq. (3), but only at a < 10^^ (!), where its excitation energy scales linearly 
with m. For higher values of a the energy remains fiat. The next two plots. 
Fig. 2(b-c), show the contributions to the sum rules as a function of logo;. As 
required by Eq. (3), the chiral soft mode saturates the sum rules in the chiral 
limit, and becomes dominant at a < 10~^. 

Figure 3 shows the results for the parameter set II. The important qualitative 
difference when compared to Fig. 1 may be seen in the bottom of Fig. (c): 
the excitation energy of the spin-isospin sound is negative now. As discussed 
in Ref. [1], the energy of a charged excitation in asymmetric medium need 
not be positive. This is because we require that the state \C) be the ground 
state for a given value of the constraint, here the isospin. A charged excitation 
connects |C) to states with different isospin, hence such an excitation may in 
fact be a deexcitations. More explicitly, we note that the u quark is lighter 
than d, because the difference of their vector self energies, p, which is nega- 
tive, overcomes the positive difference of the scalar self energies (cf. Eq (7)). 
Therefore changing m to o? is favorable energetically. Note that this negative 
excitation is still above the chemical potential for isospin, which is equal to 



5 



I^u — l^d — P + y kl + — yka + M^, hence it carries positive isospin. Thus 
the system may lower its energy by spontaneously emitting the tt^ mode. 
However, by charge conservation this has to be accompanied by the emission 
of a particle of negative charge. The strong-interaction modes (tt^, or ud pair) 
are heavy enough, so the system remains stable with respect to strong decays. 
In weak decays the accompanying particle is the electron or muon. Depending 
on the value of the chemical potential of these particles, the reaction may lead 
to a net energy gain. In that case the system is unstable with respect to weak 
interactions. 

Figures 3(c-d) show the sum rules for the parameter set II. Again, we notice 
that the spin-isospin sound mode contributes at the level of 10-20% to the 
GMOR sum rule, but with opposite sign. It dominates the isovector sum rule 
at x > 3. Figure 4 shows a qualitatively different behavior of the system close 
to the chiral hmit for the parameter set II than for the set I. We see that now 
it is the "big" 7r+ mode, not the 7r+, which assumes the role of the chiral soft 
mode. The tt^ has a finite energy in the chiral limit, and does not contribute 
to the sum rules. These are dominated by the "big" tt"*" mode at a < 10~^, 
where it scales linearly with m. 

The linear dependence of the excitation energy of the chiral soft mode on 
the current quark mass takes effect only at extremely small values of m, cor- 
responding to m^r 10 — 100 times lighter than the physical value. Hence the 
satisfaction of the laws (3) has formal rather than practical meaning. Never- 
theless, it is interesting that these laws may be fulfilled by the highly collective 
spin-isospin sound, rather than by the vacuum pion branch. In our model at 
all values of m and densities the three pionic modes saturate the current- 
algebra sum rules at the 99% level, and the contributions from the spin-isopin 
sound are big. This shows that in constructing an effective theory of mesons in 
medium one should account for such excitations, including possible colective 
modes, and that the contribution from particle-hole excitations are negligi- 
ble. It is challenging to see if similar results hold in more realistic approaches 
to nuclear matter, with nucleon degrees of freedom. In particular, the role 
of short-range correlations should be examined, since they are known affect 
strongly the excitation spectra in nuclear matter. 
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Fig. 1. Dependence of various quantities on x = ps/po for the parameter set I 
and y = 2/3: (a) Mean fields in MeV. (b) Excitation energies for the "big" modes 
(top) and the spin-isospin sound (bottom) in MeV. Also shown is the particle-hole 
production cut (shaded region), (c) The relative contributions of the three modes 
to the GMOR sum rule. The sum (dash-dot line) is indistinguishable from one. (d) 
Same as (c) for the isovector sum rule. 
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Fig. 2. Study of the chiral limit. Convention for lines as in Fig. 1(b). (a) Dependence 
of excitation energies of the three modes on the ratio of the current quark mass to 
its physical value, a. The curve for the spin-isospin sound has been multiplied by 10. 
(b) Relative contributions of the three modes to the GMOR sum rule, (c) Relative 
contributions of the three modes to the isovector sum rule. The spin-isospin sound 
becomes the chiral soft mode. 




Fig. 3. Same as Fig. 1 for the parameter set II. The spin-isospin sound mode in (b) 
has negative excitation energy. 




Fig. 4. Same as Fig. 2 for the parameter set II. The "big" mode becomes the chiral 
soft mode. 



